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Ultracold gases promise access to many-body 
quantum phenomena at convenient length and 
time scales. However, it is unclear whether the 
entropy of these gases is low enough to real¬ 
ize many phenomena relevant to condensed mat¬ 
ter physics, such as quantum magnetism. Here 
we report reliable single-shot temperature mea¬ 
surements of a degenerate ®^Rb gas by imag¬ 
ing the momentum distribution of thermalized 
magnons, which are spin excitations of the atomic 
gas. We record average temperatures as low as 
0.022(l)stat(2)sys times the Bose-Einstein conden¬ 
sation temperature, indicating an entropy per 
particle, S/N « 0.001/cs at equilibrium, that is 
well below the critical entropy for antiferromag¬ 
netic ordering of a Bose-Hubbard system. The 
magnons themselves can reduce the temperature 
of the system by absorbing energy during ther- 
malization and by enhancing evaporative cooling, 
allowing low-entropy gases to be produced within 
deep traps. 

Trapped quantum gases can be brought to impres¬ 
sively low temperatures. A single-component Bose gas 
was cooled to around 500 pK by adiabatic expansion 

[1] , and a two-component lattice-trapped Bose gas was 
cooled to 350 pK using a spin demagnetization technique 

[2] . However, the entropies achieved in these systems 
are still much higher than would be required to observe 
many-body quantum effects such as magnetic ordering 
or d-wave superconductivity of atoms in optical lattices 

[3] . Estimated critical entropies required to achieve quan¬ 
tum magnetic ordering of atoms in optical lattices include 
S/N ~ 0.3 ks and S/N ~ 0.03 ks for bosons in state- 
dependent three- and two-dimensional lattices, respec¬ 
tively [4], and similar values for Neel ordering of lattice- 
trapped Eermi gases [5], with lower entropy needed for 
less strongly interacting gases. 

In many experiments on strongly interacting atomic- 
gas systems, the low-entropy regime is reached by first 
preparing a weakly interacting bulk Bose gas at the low¬ 
est possible temperature, and then slowly transforming 
the system to become strongly interacting [6-10]. To dis¬ 
cern whether the transformation is adiabatic and to de- 
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termine indirectly the thermodynamic properties of the 
strongly interacting system, the system is returned to 
the weakly interacting regime where relations between 
temperature, entropy, and other properties are known. 
Therefore, methods to lower entropies and measure tem¬ 
peratures of weakly interacting gases are important for 
the study of both weakly and strongly interacting atomic- 
gas systems. 

In this Letter, we report cooling a Bose gas to a few 
percent of the condensation temperature, Tc, correspond¬ 
ing to an entropy per particle of 10“^ ks {ks is the Boltz¬ 
mann constant), nearly two orders of magnitude lower 
than the previous best in a dilute atomic gas [2, 11]. 
Surprisingly, we achieve this low entropy using a stan¬ 
dard technique: forced evaporation in an optical dipole 
trap, which we find remains effective in a previously un¬ 
characterized regime. The lowest temperatures we re¬ 
port are achieved at very shallow final trap depths, as 
low as 20 nK, set by stabilizing the optical intensity with 
a long-term fractional reproducibility better than 10“^. 
In addition, we demonstrate and characterize a method 
of cooling that lowers the entropy without changing the 
trap depth, possibly allowing the low entropy regime to 
be reached or maintained in systems where the trap depth 
is constrained. 

Both thermometry and cooling require a means of dis¬ 
tinguishing thermal excitations. Eor example, forced 
evaporative cooling [12, 13] depends on the ability to se¬ 
lectively expel high-energy excitations from the system. 
Similarly, thermometry of a degenerate quantum gas re¬ 
quires one to identify the excitations that distinguish a 
zero- from a non-zero-temperature gas. Both these tasks 
become difficult when S/N, or, equivalently, the fraction 
of thermal excitations Nth/N, is small [14]. Time-of- 
flight temperature measurements, in which the gas is re¬ 
leased from the trap and allowed to expand before being 
imaged, have required Nth/N of at least several percent, 
limiting such thermometry of Bose gases to T > 0.3 Tc 
and of Fermi gases to T > 0.05 Tp, where Tc and Tp are 
the Bose-Einstein condensation and Fermi temperatures, 
respectively. 

We extend thermometry to the deeply degenerate 
regime of a Bose gas by measuring the momentum dis¬ 
tribution of a small number of spin excitations, similar 
to the co-trapped impurity thermometry commonly em¬ 
ployed in Fermi gases [15-18]. Even in a highly degen¬ 
erate Bose gas with a vanishingly small non-condensed 
fraction, the minority spin population can be made di¬ 
lute enough to remain non-degenerate and thereby carry 
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FIG. 1. Magnon thermometry, (a) Magnons are created 
and decohere rapidly in a non-degenerate spinor Bose gas at 
an intermediate trap depth. Forced evaporative cooling to a 
final, variable, trap depth reduces the temperature and the 
majority gas undergoes Bose-Einstein condensation. Images 
and corresponding integrated line profiles of the momentum 
distribution of the majority (b) and magnon gas (c) are each 
shown at three different final trap depths. Line profiles are 
shown offset for clarity. A condensate obscures the momen¬ 
tum distribution of the non-condensed fraction of the majority 
gas, especially at low temperatures. In contrast, the magnons 
can have have little to no condensed fraction, allowing non- 
condensed magnons to be identified and the temperature de¬ 
termined. 


a large entropy and energy per particle. Furthermore, the 
minority spins that we use—magnons within a ferromag¬ 
netic spinor Bose-Einstein condensate—support a higher 
iVth than the majority spins because of their free-particle 
density of states [19], increasing the signal of the temper¬ 
ature measurement. Performing spin-selective measure¬ 
ments on the minority spins allows the gas temperature 
to be readily determined. 

Our procedure to measure temperatures is illustrated 
in Fig. 1(a). Experiments begin with spin-polarized ®^Rb 
in the 1 ^= 1 , tof=— 1 ) state confined in an anisotropic 
optical dipole trap and cooled to just above quantum 
degeneracy by forced evaporation to an intermediate trap 
depth. The spin quantization axis is defined by a 180 mG 
bias magnetic field. We tip the gas magnetization with a 
brief radio-frequency (RE) pulse, coherently transferring 
a small fraction (up to 15%) of the atoms primarily into 
the 11^=1,771^=0) state, creating magnons which rapidly 
decohere and thermalize. The gas is then cooled further 
by forced evaporation to a hnal trap depth where the 
temperature reaches a steady state. 

Next, we release the gas from the optical trap and im¬ 


age the momentum distribution of the magnons. Upon 
extinguishing the trap light, the gas expands rapidly in 
the most tightly confined (vertical) direction, quickly re¬ 
ducing the outward pressure within the gas. We then 
use a sequence of microwave and optical pulses to drive 
away atoms not in the |mi;’=0) Zeeman state. Finally, 
we transfer the remaining atoms to an atomic state 
{\F=2,mF=l)) suitable to two-dimensional magnetic fo¬ 
cusing [ 20 ] so that their transverse spatial distribution, 
which we image, closely reflects their initial transverse ve¬ 
locity distribution. The majority gas is probed by imag¬ 
ing |r 77 i?= 0 ) atoms immediately after the application of 
an RF pulse. 

The advantage of using incoherent spin excitations to 
measure temperature is exhibited in Fig. 1, which com¬ 
pares the momentum distribution of the majority gas (b) 
to that of the magnons (c). At T < 0.3Tc, the non- 
condensed fraction of the majority gas is obscured by the 
condensed fraction, the distribution of which is broad¬ 
ened owing to interactions and imaging resolution. In 
contrast, the temperature can be easily extracted from 
the momentum distribution of the thermalized minority- 
spin gas. Temperatures extracted from the magnons 
[Fig. 2(a)] agree with those from the majority-spin dis¬ 
tribution under conditions where both can be measured. 
At lower temperatures, the magnon thermometer allows 
measurements in a hitherto unmeasurable regime. 

Our measurements reveal that forced evaporative cool¬ 
ing produces Bose gases with extremely low temper¬ 
ature and entropy. At the lowest trap-depth setting 
(dashed circles, Fig. 2), we measure an average gas tem¬ 
perature of T = 1.04(3)stat(7)sys nK, corresponding to 
T/Tc = 0 . 022 ( 1 )( 2 ) where Tc is calculated using the mea¬ 
sured atom number, = 8.1 x 10 ®, and optical trap fre¬ 
quencies. For this gas, ksT /« 0.07, with fi the chem¬ 
ical potential, meaning that thermal excitations in the 
high density region of the condensed gas (in the major¬ 
ity spin state) are predominantly phonons. Measurement 
procedures, error estimates, and calculations of Tc and /r 
are detailed in Methods. 

We observe that evaporative cooling is highly efficient 
with respect to particle loss in the regime ^/ks < T < Tc 
[Fig. 2(b)]. When T < /i/Zcs, the cooling efficiency re¬ 
duces, consistent with the fact that the chemical poten¬ 
tial accounts for a non-negligible amount of the energy 
carried away by each atom lost to evaporation. The rele¬ 
vance of the chemical potential in evaporation is also indi¬ 
cated in Fig. 2(a) by the observation that 77 = [//(fc^T), 
the ratio of trap depth to thermal energy, increases at 
lower trap depths, and that the temperature at a fixed 
trap depth depends on the number of atoms. In an evap- 
oratively cooled gas, the temperature responds to the 
effective trap depth, C/es = U — the potential energy 
depth minus the chemical potential, rather than the trap 
depth alone, and in the regime T < 11 /ks, the difference 
U — UeS = /7 is significant. 

We calculate the Bogoliubov energy spectrum of the 
confined quantum degenerate gas, including the effects 
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FIG. 2. Thermometry results. Two runs (blue circles, green squares) differ in initial atom number, (a) T and (b) T/Tc are 
measured at various optical trap depths. Lower T and T/Tc are achieved in runs with larger initial atom number. In (a), ther¬ 
mometry using the majority spin population (light gray diamonds) agrees with thermalized magnon thermometry extrapolated 
to the zero-magnon values (circles and squares). Error bars show statistical uncertainty of extrapolated temperature (vertical) 
and systematic uncertainty in trap depth (horizontal). Thin diagonal gray lines show contours of rj, the ratio of trapping poten¬ 
tial depth to temperature. Thick gray lines show calculated Tc and ^/ks in (a) and fi/ksTc, calculated using the frequencies 
of the fes X 60 nK deep trap, in (b). (b) Efficiency of evaporation only (circles, squares), with error bars showing statistical 
uncertainty of the extrapolated T/T^, is compared to magnon-assisted evaporation (“x”). Solid lines connecting points are 
guides to the eye, with a steeper slope indicating more efficient evaporation. Dotted lines connect corresponding points at the 
same trap depth. Inset: the magnon-free temperature is determined by extrapolating single-shot temperatures (“x”) vs. total 
number of magnons imaged, here at trap depth of 800 nK, to the zero-magnon limit (circle). A similar extrapolation determines 
the value of T /Tc achieved by evaporation only. 


of trap anharmonicity, and find its entropy at equilib¬ 
rium (see Methods) to be S/N = 1 x 10“^ ks, the low¬ 
est value ever reported for an atomic gas. By com¬ 
parison, the 500 pK Bose gas reported in Ref. [1] has 
S/N = 3.6fcB(T/Tc)^ ^ 1.5 ks using relations for a non¬ 
interacting gas at T/Tc ^ 0.75. Other reported values 
include S/N ^ 0.05 fcs at the center of a resonantly in¬ 
teracting Fermi gas [11], S/N = 0.27fcB[21] or ^0.1 A:b[2] 
for bosons in a lattice, and T/Tc 0.15 in a double¬ 
well Bose-Einstein condensate [22]. Thermometry based 
on imaging incoherent phonons indicated a temperature 
around 7 times higher than reported here in a ®^Rb con¬ 
densate of similar density [23]. 

Having magnons present during evaporative cooling re¬ 
duces the temperature of the trapped gas by increas¬ 
ing the efficacy of evaporative cooling. Forced evapo¬ 
rative cooling from a trap with an effective trap depth 
UeS 3> ksT has a cooling power proportional to the 
number of thermal excitations with excitation energies 
above C/eff- In a weakly interacting single-component de¬ 
generate Bose gas, the number of thermal excitations is 
determined by the temperature and is independent of the 
total particle number, fixing the evaporative cooling rate. 
By seeding the gas with additional spin excitations at 
constant total particle number, the total number of ther¬ 
mal excitations, and thus the evaporative cooling power, 
increases. 


We observe that, for T > /z//cb, this magnon-assisted 
evaporation leads to T /Tc lower than that reached by 
single-component evaporation at the same final particle 
number (“x” points in Fig. 2(b)). When T < t/^b, we 
do not find magnons to reduce the temperature, perhaps 
because of the small number of non-condensed magnons 
and the decrease in evaporative cooling efficiency in gen¬ 
eral. Magnon-assisted evaporation is also evidenced by 
the fact that the per-atom loss rate is higher for the 
magnon population than for the spin-majority popula¬ 
tion: For both populations, the atom loss reflects the 
selective evaporative loss of thermal excitations. For the 
spin-majority atoms, thermal excitations are just a small 
fraction of the total population, whereas they represent 
a much higher fraction of the magnons. 

Unbiased thermometry of the majority gas is per¬ 
formed by extrapolating measurements with various 
numbers of magnons to the zero-magnon limit (Fig. 2(a) 
inset). In addition, by examining evaporative cooling in a 
state-dependent optical trap, we have confirmed that the 
temperature indicated by the magnons varies with the 
trap depth of the spin-majority gas, even for a fixed trap 
depth of atoms in the |mp’=0) state. This observation 
gives further confidence that the magnon thermometer 
accurately records the temperature of the majority gas. 

The same characteristics that make dilute magnons a 
good thermometer—high relative energy and entropy per 
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FIG. 3. Cycled decoherence cooling, (a) Magnons are created at the final trap depth and cool the gas as they thermalize. 
Thermalized magnons can be removed from the trap and the cooling process repeated before measuring the temperature by 
imaging the momentum distribution of the thermalized magnons. (b) Each non-condensed magnon removed from the trap 
takes energy from the gas. Decoherence cooling trajectories are plotted versus cumulative number of magnons removed. Closed 
circles show temperatures after a single cycle of decoherence, but with varying numbers of magnons. Some additional magnon- 
assisted evaporative cooling may also be present. Open triangles show repeated cycles of magnon creation, thermalization, and 
purge, with representative error bars on first and last triangles giving statistical error over several repetitions. Solid lines show 
zero-free-parameter theory predictions assuming each non-condensed magnon removes 3 ksT energy. Solid patches indicate 
the range of predictions included within the uncertainty in the cumulative number of magnons removed (we image 75% ± 10% 
of the magnons present). Inset: the amount of energy removed by each magnon, in units of ksT, is plotted versus the ratio 
of n/ksT, where each point is based on an estimate of the slope over 4 consecutive triangle points. We observe that the net 
energy carried away by impurities vanishes when T < ^i/ks- (c) Efficiency of evaporation only (black circles) is compared to 
cycled decoherence cooling (open triangles), with error bars representing statistical uncertainty. The solid line between circles 
is a guide to the eye. Dotted lines connect corresponding points at the same trap depth. Compared to evaporation, decoherence 
cooling can reach considerably lower T/Tc at the same trap depth in the regime T > fi/ksT and T/Tc < 0.9. 


particle—enable them to cool without lowering the trap 
depth, beyond what can be achieved by magnon-assisted 
evaporation. Spin excitations created in a degenerate 
Bose gas decrease its temperature by a process known as 
decoherence [24] or demagnetization [25, 26] cooling. Im¬ 
mediately after the RF pulse is applied, the magnon pop¬ 
ulation has the same energy and momentum distribution 
as the initially polarized degenerate Bose gas, with a large 
condensed fraction, which carries no energy above the 
chemical potential, and a small normal fraction, which 
carries on the order of ksT of energy per particle. Upon 
thermalization, the normal fraction within the magnon 
gas increases, bounded from above by its critical num¬ 
ber for condensation, and the energy and entropy of the 
magnons increases. The energy gained by the magnons 
is energy lost by the majority gas, as the entire process, 
apart from the slight effects of magnetic field inhomo¬ 
geneities and dipolar interactions, occurs at constant en¬ 
ergy [27]. 

While the overall entropy rises during thermaliza¬ 
tion, the entropy per particle of the majority gas may 
decrease. Considering a non-interacting harmonically 
trapped gas, we calculate that decoherence cooling would 
reduce T/T^ of the majority gas only when T/T^ < 
^ (3C(3))/(4^(4)) = 0.94, consistent with our observa¬ 


tions. 

To demonstrate the capabilities of decoherence cooling, 
we introduce a variable number magnons at the final trap 
depth and allow them to evolve toward thermal equilib¬ 
rium with the assistance of a weak magnetic field gradient 
of 0.2-2 mG/cm (for more on magnon thermalization, 
see Supplemental Discussion). Then, as illustrated in 
Fig. 3(a), we either measure the temperature or expel the 
magnons using a spin-selective process, leaving a colder 
spin-polarized condensed gas. The decoherence cooling 
process can then be repeated to further cool the sample, 
with the number of magnons, magnetic field gradient, 
and thermalization time optimized at each repetition. A 
single pass of decoherence cooling, using pseudo-spins in 
a degenerate ®^Rb gas, has previously been shown to 
slightly reduce the temperature, but not T/Tc, of the 
resulting mixture [24]. 

Whereas the thermal energy kgT achieved by forced 
evaporative cooling of a single-component gas is limited 
to around 1/t] ~ 1/10 of the trapping potential depth, the 
temperatures reached by cycled decoherence cooling are 
far lower. For example, we produced a 20 nK gas within a 
fcs X 600 nK deep potential, corresponding tol /?7 = l/30, 
and reduced T/Tc substantially, from 0.58 to 0.37 [cyan 
triangles. Fig. 3(b) and (c)]. On the other hand, de- 
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coherence cooling is generally less efficient in terms of 
atom loss. An atom lost through evaporative cooling re¬ 
moves an energy ~ vyfcsT, whereas a magnon brought to 
equilibrium and forcibly ejected only removes 3 ksT for 
a non-interacting gas, consistent with our observations 
when ksT 3> At lower temperatures, we expect each 
impurity to remove less energy, 1 to l-SfesT, as interac¬ 
tions and quantization of the trap modes become more 
important. Even so, cooling would be expected, in the 
ideal case, to become increasingly effective as the heat 
capacity drops. In contrast, we hnd decoherence cool¬ 
ing to be ineffective at reducing both T and T/Tc in the 
low temperature regime ksT < /x (Fig. 3(b, c), inset and 
red triangles), perhaps owing to energy deposited in the 
trapped gas during the magnon expulsion and by heating 
processes, which become more significant as thermaliza- 
tion times increase. 

In principle, the process of creating the spin mix¬ 
ture employed in decoherence cooling can heat the gas 
by changing its interaction energy. However, neglecting 
small magnetic dipolar interactions, the magnon is a gap¬ 
less Goldstone excitation corresponding to the rotational 
symmetry of real spins [19]. Accordingly, the s-wave scat¬ 
tering lengths that characterize the interactions between 
two majority-spin atoms, and between a majority 

and minority-spin atom, a_i^ 0 ; are identical. In contrast. 


spin excitations involving a different hyperfine state can 
deposit energy into the condensate due to a mismatch of 
s-wave scattering lengths. 

In this work, we have demonstrated the use of spin ex¬ 
citations within a highly degenerate Bose gas to measure 
and reduce the temperature. With our improved ther¬ 
mometry, we observe a record low entropy of 1 x 10“^ ks 
per particle for a spin-polarized Bose-Einstein conden¬ 
sate, illustrating the power of standard cooling tech¬ 
niques to reach the long-sought low entropy regime. 
Magnon-assisted evaporation and cycled decoherence 
cooling should work for other spin or pseudospin mixtures 
of Bose gases so long as the gas thermalizes at constant 
minority particle number, providing a flexible platform 
for creating or maintaining low entropies in deep traps. 
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METHODS 

A. Optical trapping 

Experiments were performed on ®^Rb gases in a single¬ 
beam optical dipole trap [Ml]. Evaporative cooling was 
realized by gradually lowering the depth of the opti¬ 
cal trap. The optical trap was formed by light with 
a wavelength of 1064 nm, brought to a cylindrical fo¬ 
cus at the location of the atoms. The optical trap was 
highly anisotropic, with trapping frequencies having a 
typical ratio 1:15:140, with the tightest confinement in 
the vertical direction. At low optical powers, the optical 
dipole trap is influenced by gravity so that the ratio of 
trap frequencies differs somewhat. Trap frequencies were 
measured empirically at several optical powers, extrapo¬ 
lated between measurements, and confirmed by compar¬ 
ing with a simple model of the trapping potential that 
accounts for the Gaussian focus and force of gravity. 

The bulk of our experiments were performed with lin¬ 
early polarized trap light, which imparts equal dipole 


force on all Zeeman sublevels of F=1 ®^Rb . To ver¬ 
ify that our temperatures vary with the trap depth of 
the majority gas, rather than that of the minority spins, 
we also performed experiments with circularly polarized 
trap light, which applies a differential force to atoms in 
different Zeeman sublevels. Using this circularly polar¬ 
ized trap, we vary the trap depth of the majority gas 
without affecting that of the |mi?=0) atoms by preparing 
the majority gas in either of the \mp=+l) or \mp=—l) 
states. 

The lifetime and heating rate of the atoms in the op¬ 
tical trap are consistent with spontaneous scattering of 
the far-detuned trapping light. 


B. State purification and preparation 

About one millisecond after extinguishing the opti¬ 
cal trapping light, we initiate the imaging sequence by 
purging atoms in the \F=l,'mp=—l) state—the major¬ 
ity component of our gas in most of our experiments— 
using a series of alternating microwave pulses (resonant 
with the \F=l,mp=—\) —> \F=2,mp={)) transition) 
and light pulses (resonant with the F=2 hyperfine man¬ 
ifold). 

The number and strength of the optical/microwave 
pulses is optimized to reliably expel all of the 
1^=1, mp=—l) atoms with minimal impact on the num¬ 
ber of atoms in the \F=l,mp=0) state. Regardless, 
the purge process removed between 10% and 40% of 
the |F=1 ,tof= 0) atoms, depending on the number of 
\F=l,mp=—l) atoms being purged and other experi¬ 
mental parameters, but it does not visibly impact the 
momentum of the atoms that remain. We estimate the 
number of magnons from the images by accounting for 
this loss, allowing considerable uncertainty. 

Finally, we employ a calibrated microwave sweep to 
transfer atoms in the \F=l,mp=0) state to the mag¬ 
netically trappable |F=2, mp=l) state with greater than 
95% efficiency. 


C. Magnetic focusing 

The magnetic focusing lens takes the form of a mag¬ 
netic potential with negligible curvature along the verti¬ 
cal (imaging) axis and weak, harmonic curvature in the 
image plane. The resulting magnetic trap causes atoms 
with momenta p in the image plane to converge on the 
real-space point Xp = plmio after a time tfoc = 27r/4a;, 
where ra is the atomic mass and w the trap frequency 
for atoms in the \F=2,mp=l) state. The vertical gra¬ 
dient is selected to cancel the effect of gravity, and the 
trap frequency oj is chosen such that the thermal atomic 
momentum-space cloud is resolvable by our imaging sys¬ 
tem at the lowest accessible temperatures while still al¬ 
lowing us to image temperatures near Tc. 
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Uncertainty in the parameter oj is the largest source 
of systematic error in the temperatures reported in this 
work. Our measurements found w = 27rx2.88±0.09 rad/s 
and the relation Xp = p/mu! was verified by comparing 
the amplitude and phase of center-of-mass oscillations of 
our gas imaged both in situ and after application of the 
magnetic focusing lens. 

Time-of-flight expansion of the gas in the unconfined 
vertical direction along with misalignment of the trap, 
imaging, and focusing axes can lead to several types of 
aberration. For example, any projection of the vertical 
extent of the gas into the imaging plane will manifest 
as additional apparent momenta in the imaged column 
density. Such aberrations can lead to an overestimate of 
the temperature of the gas and are most pertinent at low 
temperatures. Trap, imaging, and focusing axes are all 
aligned to gravity in our system, and can be verified by 
imaging highly degenerate gasses, which expand negligi¬ 
bly in the direction of weakest optical confinement. 


D. The imaging system 

After magnetic focusing, the atoms are imaged with a 
40 to 100 ps pulse of light resonant with the F=2 —)■ F=3 
cycling transition. The magnification of our imaging sys¬ 
tem was calibrated using an optical micrometer target 
that was placed in a plane equivalent to the plane of the 
atoms, with respect to the imaging system. This cali¬ 
bration was repeated several times with fractional uncer¬ 
tainty of 1% and is included in the systematic error of 
our thermometer. 

The resolution of our imaging system, approximately 
8 /rm, corresponds to a temperature of about 0.2 nK. We 
do not correct extracted temperatures for finite imaging 
resolution. 


E. Details of temperature fitting procedure 

Atom column momentum densities Up are fit to a Bose- 
enhanced momentum distribution Up = b + Agj{e^), with 
gj the Bose function (polylogarithm) of order j and the 
argument 

^ (yPw Pwo'j {pt PtO^ /^bTi. 

For non-degenerate gasses, a = p/ksT is a free parame¬ 
ter of the fit, along with the zero of momentum {pwO,Pto)j 
the background level b, peak level A, and temperatures 
Tu, and Tt along the weak and tight in-plane axes of the 
optical trap, respectively. For degenerate gasses, 0 = 0. 

Generally, T-u, = Tt within the error of the fit, however 
at low temperatures, several effects cause Tt be an un¬ 
reliable estimate of the gas temperature. Although the 
condensate expands primarily along the unfocused verti¬ 
cal direction, at the very low temperatures reported in 
this work, condensate expansion along the more tightly 


confined in-plane direction can be manifest. Also, the 
tilt of the trap in the tight direction is more difhcult 
to calibrate, leading to a (fictitious) systematic upward 
shift in the apparent Tt by the mechanism explained in 
Sec. C. Finally, at the lowest temperatures reported in 
this work, the semiclassical condition ksT ^ bio does 
not hold along the tight axis as it does along the weak 
axis {h{oJw,oot}/kB ~ {0.05, l}nK). Thus, in this work, 
we estimate temperatures by T = T^, alone. 

All fits exclude regions of the column density that in¬ 
clude condensed atoms, and each fit was performed many 
times, varying both the size of the excluded central region 
and the order j € [0.5, 2] of the Bose function in order 
to look for systematic shifts in the temperature. The 
high-momentum tails of the momentum distribution Up 
are insensitive to the particular Bose function employed 
as —)■ 0). We used only hts where the 

variation of the temperature with respect to the size of 
the exclusion region and the order of the Bose function 
was negligible compared to other sources of systematic 
error. 

The uncertainty with which individual fits estimate the 
temperature is consistent with fundamental sources of 
noise, photon and atom shot noise, and the short-time 
(same day) shot-to-shot variation of estimated tempera¬ 
tures was consistent with the single-shot noise-limited er¬ 
ror estimated by the fitting routine. This justifies the use 
of an uncertainty-weighted average of many independent 
measurements to compute the precise low temperatures 
reported in this work. 


F. Calculation of the entropy per particle 

Calculations of the entropy per particle cannot rely 
on the local density approximation because the thermal 
phonon wavelength is far larger than the transverse ex¬ 
tent of the condensed gas. We can, however, apply a 
one-dimensional local density approximation along the 
long axis of the condensate. 

To calculate the entropy, we numerically solve for the 
Gross-Pitaevskii ground state and Bogoliubov spec¬ 
trum {ci} of small amplitude excitations in a weakly- 
interacting ®^Rb Bose gas with chemical potential p in 
a trap with confinement in 2 dimensions, including the 
effects of gravity and anharmonicity [M2], and no con¬ 
finement in the 3rd dimension. To do this calculation, 
we discretize the 2-dimensional Gross-Pitaevskii equa¬ 
tion and Bogoliubov equations on a Lagrange mesh based 
on the Hermite polynomials [M3, M4]. The Bogoliubov 
equations include the substitution p ^ p — h'^k^/2m, 
for relevant values of the wave vector k, to account for 
excitations in the unconfined longitudinal direction. We 
verify that our calculation of the excitations {ci} are ac¬ 
curate by comparing to analytic results in regimes where 
they are available, for example in the harmonic non¬ 
interacting limit and in the Thomas-Fermi limit with 
cylindrical symmetry. 



The entropy and number of atoms per unit length can 
then be calculated using the equilibrium relations 


S 

L 

N 

~L 


(Si — m) /^bT 


log (l - 


A^c + E 

I 


1 


G. Calculating Tc and p 

The estimates of Tc and fj, reported use standard for¬ 
mulas for Bose gases in the non-interacting and Thomas- 
Fermi limits, respectively [M5]. The effects of trap an- 
harmonicity and finite interactions on both fx and were 
considered and are on the order of the systematic uncer¬ 
tainty in our temperature measurements. 


with Nc = 'I'*11'. By fixing T and varying /i (dl and {ci} 
are recalculated for each value of fx), we can consider the 
entropy and number per unit length as functions of ix: 
f (m) and 

We use the local density approximation and Thomas- 
Fermi approximation along the long axis of our trap (with 
trap frequency Wz) to write the total entropy of the gas 
as 


5 = 



S 

L 



dz, 


and the total number of atoms as 


N 



^Mmax 


1 2 2 
2 ^ 


dz, 


where /imax is chosen to yield the experimentally mea¬ 
sured value of N. The entropy per particle of the gas is 
then S/N. 

At the lowest trap depths, the rate at which the 
magnons are lost by forced evaporation becomes com¬ 
parable to the rate at which non-condensed magnons are 
generated through thermalization. We speculate that the 
number of thermal excitations in the spin-majority gas 
under these conditions is similarly lower than at the ob¬ 
served temperature in equilibrium. As such, our estimate 
of S/N = 1 X 10~^kB may be an underestimate of the 
true steady-state entropy. 


H. Theoretical limits for decoherence cooling 


Considering an ideal Bose gas in a harmonic trap with 
geometric-mean frequency w, the following relations for 
the saturated thermal atom number A^th: energy E, and 
heat capacity C, at the temperature T, apply: 


A^th = C(3) 
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Suppose we transfer dN > 0 atoms to the minority 
spin state (where dN is far below the critical number for 
magnon condensation), let them thermalize, and eject 
them. To lowest order, the majority gas loses energy 
dE = —SkBTdN, and thus the temperature changes by 
dT = dE/C = —3{kBT)/CdN. The non-condensed 
fraction changes as, using Eqs. Ml, 
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Thus, we must have Nth/N < (3^(3)) / (4^(4)) = 0.83 to 
reduce the non-condensed fraction through such cooling. 
This corresponds to, for the case under consideration, 
T/Tc < (0.83)1/3 = 0.94. 
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DISCUSSION 
Magnon thermalization 

We found that applying a magnetic field gradient was 
necessary to achieve thermalization between initially co¬ 
herent, quantum degenerate spin components, but can 
be chosen to be small enough to impart negligible energy. 
A gradient of 0.2 mG/cm, typical in our experiments on 
decoherence cooling, corresponds to an energy of approx¬ 
imately h X 10 Hz = fcs X 0.5 nK across the longest dimen¬ 
sion of our gas and was sufficient to produce a saturated 
non-condensed fraction of the \mF=0) state in 100-2000 
ms, depending on the temperature of the gas. Thermal¬ 
ization of magnons created out of a degenerate gas in two 
representative cases of ksT/fi < 1 and ksT/pi > 1 are 
shown in Fig. SI (a) and (b). Fig. SI (c) shows the ex¬ 
tracted temperature of magnons created at a higher trap 
depth (in a non-degenerate gas) reaching a steady state 
following the final phase of evaporative cooling, for a rep¬ 
resentative case ksT/fj, < 1. Temperatures extracted 
from magnons created from a non-degenerate gas before 
completing evaporation and from magnons that are cre¬ 
ated in the degenerate gas after evaporation agree with 
each other, when the same experimental parameters are 
employed, and when extrapolated to the zero-magnon 
limit. 

The physics behind the decoherence and thermaliza¬ 
tion of spin excitations, and the variables that affect it, 
remains a compelling direction of future research. Pre¬ 
vious experiments on degenerate F=1 ®^Rb gases have 


found magnetic field inhomogeneity to facilitate deco¬ 
herence and thermalization of spin populations, typically 
on timescales of several lOs-lOOs of milliseconds [S1-S3]. 
We have previously observed coherence between magnons 
in different momentum states out to several hundred 
milliseconds [S4] with a negligible gradient. However, 
none of these experiments claim to have operated in the 
highly-degenerate phonon-dominated regime ksT < fj.. 
With very small gradients, we find the zero-momentum 
magnons created in this work can remain coherent for 
many seconds. In strongly-interacting Fermi gases, de¬ 
coherence of spin excitations has been well characterized 
and is limited by a fundamental lower-bound on the spin 
diffusivity, h/m [S5, S6]. In this context, a gradient is 
essential to drive the diffusive spin transport. 

When ksT < pt, the thermal excitations of a ho¬ 
mogeneous Bose-Einstein condensate are phonon-like, 
and thus thermalization of the free-particle-like magnons 
with the majority gas requires collisions between quasi¬ 
particles with very different dispersion. Such interactions 
are expected in the first order beyond the Bogoliubov 
mean-field theory (the Beliaev theory) [S7], though to 
our knowledge no one has produced an estimate or direct 
measurement of the thermal phonon-magnon cross sec¬ 
tion. We speculate that thermalization of the magnons 
occurs primarily at the edges of the condensate where 
low condensate density implies a better match between 
majority spin and minority spin dispersion. In addition 
to facilitating decoherence, the applied gradient may aid 
in thermalization by driving the transport of low-energy 
spin excitations to the edge of the condensate. 


[51] Chang, M.-S., Qin, Q., Zhang, W., You, L. & Chapman, 
M. S. Coherent spinor dynamics in a spin-1 Bose con¬ 
densate. Nat. Phys. 1, 111-116 (2005). 

[52] Lewandowski, H. J., McGuirk, J. M., Harber, D. M. & 
Cornell, E. A. Decoherence-Driven Cooling of a Degener¬ 
ate Spinor Bose Gas. Phys. Rev. Lett. 91, 240404 (2003). 

[53] Kronjager, J. et al. Evolution of a spinor condensate: 
Coherent dynamics, dephasing, and revivals. Phys. Rev. 
A 72 , 063619 (2005). 

[54] Marti, G. E. et al. Goherent Magnon Optics in a Ferro- 


* ryan@efrus.com 

1 Present address: JILA, National Institute of Standards and 
Technology and University of Colorado, Boulder, Colorado 


magnetic Spinor Bose-Einstein Condensate. Phys. Rev. 
Lett. 113, 155302 (2014). 

[S5[ Bardon, A. B. et al. Transverse demagnetization dynam¬ 
ics of a unitary Fermi gas. Science 344, 722-4 (2014). 
[S6[ Koschorreck, M., Pertot, D., Vogt, E. & Kohl, M. Univer¬ 
sal spin dynamics in two-dimensional Fermi gases. Nat. 
Phys. 9, 405-409 (2013). 

[S7[ Phuc, N. T., Kawaguchi, Y. & Ueda, M. Beliaev theory 
of spinor BoseEinstein condensates. Ann. Phys. (N. Y). 
328, 158-219 (2013). 


80309-0440, USA; edward.marti@jila.colorado.edu 



2 





Thermalization time (s) 


Thermalization time (s) 


c 



Time after evaporation ramp (s) 


FIG. SI. Magnon thermalization. Magnons created in a degenerate gas with (a) ksT ~ 2^ and (b) ksT ~ /i/3 thermalize at 
trap depths of approximately fcs x 1,200 and 35 nK, respectively, in the presence of a weak magnetic field gradient. Throughout, 
traces (colored points and connecting lines) differ in initial number of magnons. In (a), the momentum distribution of the non- 
condensed magnons can be determined as soon as they are created. The indicated temperature reaches a steady state and 
the number of non-condensed magnons saturates after a few hundred milliseconds. The steady-state temperature depends on 
the number of magnons created owing to the effects of decoherence-driven cooling. In (b), temperatures and the number of 
non-condensed magnons cannot be determined until after nearly one second of thermalization, at which time the temperature 
has already reached steady state, (c) Magnons created before the final forced evaporation ramp are incoherent throughout the 
whole ramp and reach a steady-state temperature ksT ~ /r/3 after two to three seconds at the final trap depth of approximately 
ks X 25 nK. 









